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The generalized Kelvin-Helmholtz (KH) and Miles mechanisms of the water wave generation by 
wind are investigated for two-layer piece-wise linear model of the wind profile. It is shown by 
asymptotic expansions in small air-to-water density ratio that two mechanisms of the instability 
operate in quite different scales. Miles' short waves are generated by weak winds, in particular. 
Miles' regime is responsible for initiation of the instability at the minimum wind speed, while the 
generalized KH regime dominates at strong winds and raises moderately-short waves. 
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I. INTRODUCTION 

The current study is motivated by the attempt to un- 
derstand and explain qualitatively water wave generation 
by winds. In particular, it aims at shedding light on 
and gaining a deeper understanding of the relative roles 
played by the Kelvin-Helmholtz (KH) and Miles mecha- 
nisms in various regimes of wavelength and wind-strength 
values. 

The Kelvin-Helmholtz and Miles theories are tradition- 
ally two principal models which describe the energy and 
momentum transfer from wind to waves. The classical 
KH instability (KHI) theory considers a tangential dis- 
continuity between uniform flows of air and water [l|. 
Exact solutions of the linear KHI can be obtained due 
to the constant profiles of the unperturbed velocities of 
both air and water. Consequently, the growth rate of the 
KH mode is evaluated explicitly. It predicts a minimum 
wind speed necessary to make the waves grow against 
the stabilizing effects of gravity and surface tension. The 
quasi-laminar Miles model takes turbulence effects into 
account through the mean profile of the unperturbed 
wind which continuously varies with the height above the 
air- water interface [3] . The characteristic feature of the 
Miles model is the concept of resonant wind- wave interac- 
tion, which assumes an intensive exchange of energy and 
momentum between the collective wave motion and the 
air particles within the critical layer (a narrow vicinity 
of the point within the air flow where the perturbation 



phase velocity is equal to the wind speed). 

In the limit of small air-to-water density ratio e^ — 
Pa/pw, Miles' model yields the following threshold val- 
ues of the minimum wind speed and corresponding wave- 



length for the onset of surface waves 
0.23m/s, A(^'') = 2nLc ~ 1.7 • 10" 
is the capillary length; cr, g, Pa 



u): 



(cr) 



and Pu 



n/ML 



= \Jol{pwg) 
cr, 5, Pa ana Pu, are the air- 
water surface tension, gravity acceleration, air and wa- 
ter density, respectively; a j p^ « 7.2 • lO^^m'^/s^, L^, ~ 
2.7 • 10~'^r?T,). In the leading order in e the correspond- 



ing values for the KHI are Ua 
A(='-) = 27rLc ~ 1.7 • lO-^TO. 
critical wave lengths (A'^'^''^ = 



(cr) 



^\/2gLc ~ 7m/ s, 
The same value of the 
2ttLc) in both models at 

{cr) 
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quite different threshold velocities {Ua ~ e^V^ff^c and 
Ua ^ e^^y/2gLc) sometimes leads to the misconcep- 
tion of the Miles and KH mechanisms of surface wave 
instability as alternative and competing mechanisms of 
wave generation in the same wave-length range. Thus, 
according to a widespread opinion, the Miles model is 
an improvement of the KH model (e.g. Q), and the KH 
regime has no physical meaning in the context of the sur- 
face wave generation. Such a point of view may be traced 
to 01 where it is written that 'The model to be developed 
here improves on the Kelvin-Helmholtz model by allow- 
ing for distributed, rather than concentrated vorticity...', 
in spite of the quite definite point of view later formulated 
in Q: '...the KH mechanism of instability still should be 
physically significant, albeit not responsible for the initial 
formation of water waves.' An expression of such a belief 
may be found in 5] who have artificially superposed a 
velocity jump at the water-air interface on a continuous 
wind profile. That has been done in order to allow the ex- 
citation of the KHI at short wavelengths range similar to 
that of the Miles mechanism, namely, at the order of the 
capillary length. In fact, the unphysical negative phase 
velocity of the KH mode, exhibited in [5] for a smooth 



wind profile, is a direct result of such artificially intro- 
duced concentrated vorticity. As is shown in the current 
work, the KHI does occur for continuous profile but for 
longer wavelengths and higher wind speeds than in Miles' 
regime, and the different scaling of both the characteris- 
tic wavelength and wind speed with e provides an impor- 
tant criterion to distinguish between the generalized KH 
and Miles mechanisms. Note that for wave lengths much 
larger than the capillary length, the classical KH model 
yields for the threshold wind speed Ua' « e~^^/gjk, a 
value that depends on the wave length {k — 27r/A). This 
leads to the proper characteristic wind speed for the KH- 
like regime Ua ~ e'^^^VgL^ which is achieved at wave 
lengths of the order of A ~ e^^Lc (see below Section 4, 
where it was established that accounting for the surface 
tension effect for prediction of the minimum wind speed, 
which is correct for Miles' regime, is asymptotically small 
at KH scales of wave lengths, and consequently has no 
physical meaning in the classical KH model. 

Although various aspects of the role of the KHI in wa- 
ter wave generation have been discussed previously [J]- 
[9|, and conditions necessary for the dominant role of 
the KHI have been well recognized {Ua ^ \/9La with 
X'^ La, where La is the characteristic length of the wind 
profile inhomogeneity, the roughness length, which is typ- 
ically of ~ Lc) [7], the dominant role of the Miles mech- 
anism is sometimes expanded to relatively long waves 
(A ^ La) and strong winds {Ua ^ y/gLa). In fact, as 
is shown in the current work the KH and Miles mech- 
anisms act in quite different scales of both characteris- 
tic wind speeds and wave lengths (see Sections 3 and 
4, where the corresponding orders in e have been estab- 
lished) . One should expect the KH model with a concen- 
trated vorticity to be more satisfactory for strong wind 
speeds and perturbations with wavelengths that are much 
larger than the roughness length of the smooth profiles. 
In contrast, the Miles model is more adequate for predic- 
tion of stability characteristics for low wind speeds and 
wave lengths. Finally, in addition to the large difference 
in the scales of the perturbation wavelength and wind 
velocity for the two modes, the growth rate of the KH 
mode is much larger than that for Miles' mode (see be- 
low Sections 3 and 4). 

Although the numerical solution of the problem that 
describes the wave generation by a wind of any smooth 
profile may be easily developed for a wide range of param- 
eters (see e.g. @), the simple piece- wise linear (PWL) 
models with low number of the layers N are widely used 
in geophysics and astrophysics 10] -|14|. In particular, 
the two-layer PWL approximation for the wind speed 
profile is widely employed in order to simplify the study 
of the wind instability problem and simultaneously to 
conserve the main features of the wind-shear effects (see 
e.g. 0, il|-[3). The advantage of the PWL approx- 



imations with low TV is in the explicit solutions for the 
eigenfunctions which results in a polynomial dispersion 
relation between perturbation frequency and wave num- 
ber, and allow to gain deeper understanding of the mech- 



anisms that govern the instability, and the corresponding 
principal scales for wavelengths and wind speeds, mag- 
netic fields etc. In general, such a description should be 
considered as the simplest schematic modeling of the dif- 
ferent shear effects rather than an approximation to any 
smooth wind profile that may be justified by the uncer- 
tainty in our knowledge of the wind profile near the waved 
interface. The validity of two-layer PWL approximation 
of the smooth wind profiles will be discussed in the last 
section of the present study, while the convergence of the 
results of the general PWL (Rayleigh) method with in- 
creasing N to those of the corresponding smooth profiles 
of the wind has been established in [17|-[18|. 

Simultaneous study of the the Miles and KH mecha- 
nisms of instability is complicated due to the necessity 
to solve the problem in two quite different regimes (with 
different wavelength of the perturbations and the unper- 
turbed wind speed). For instance, both 6] and [16] apply 
the same two-layer PWL approximation for the wind pro- 
file to different regimes. Fabrikant and Stepanayantz |la | 
investigate the resonant Miles mechanism of instability, 
which is characterized by a near resonance distribution 
of the growth-rate coefficient localized in the vicinity of 
two (for cr > 0) resonant wave numbers. The model by 
Fabrikant and Stepanayantz in the Miles' region of pa- 
rameters (low wave length and wind velocity) yields nu- 
merically the critical parameters for instability that are 
in fair agreement with Miles' result: Ua ~ 0.27to/s for 
wavelength A^'^''-' w 1.7 • 10~^?ti, while the KH mode was 
out of their consideration. However, the growth rate co- 
efficient of the Miles mode in the two-layer PWL approx- 
imation is significantly larger than that for the smooth 
profiles in [5]. As distinct from J16j|, Capponi et al. Q 
have solved the problem in the KH parameter regime 
that has no-resonance nature. Their model yields that 
the instability starts from vanishing wind velocity and 
wave length (in the scales of the KH mode of relatively 
long waves and strong winds). Two marginal curves ap- 
pear instead of one in the classical KH model, when the 
order n of the dispersion relation increases from n = 2 in 
the classical piece-wise constant KH approximation (vor- 
tex sheet instability) to n = 3 in the two-layer PWL 
KH-like approximation. One of the marginal stability 
curves corresponds to a generalized KH mode modified 
by the wind-shear effects, while the other marginal sta- 
bility curve is completely induced by wind-shear effects 
on the KH-like mode. That shear-induced instability has 
been erroneously interpreted in [H as the Miles mode (see 
Section 4 of the present study). 

The present study of the surface wave instability is car- 
ried out by asymptotic expansions in small air-to-water 
density ratio and applied to two-layer PWL wind ap- 
proximation. The application of the asymptotic expan- 
sions to the simplified PWL model for the wind profile 
leads to explicit solutions for both the KH regime which 
takes into account shear effects and Miles' regime. The 
present study demonstrates that both mechanisms are re- 
sponsible for wind instability, but the Miles and KH-like 



regimes are distinguished by the scales in the wavelengths 
and wind speeds. 

The paper is organized as follows. The dispersion 
relation for two-layer PWL approximation of the wind 
speed profile is presented in the next Section. The Miles 
and KH-like mechanisms of wave generation are analyzed 
asymptotically in the small density ratio in Sections 3 and 
4, respectively. Summary and Discussion are presented 
in Section 5. Asymptotic relations for the Miles regime 
are derived in Appendix A. 



II. THE DISPERSION RELATION FOR 
TWO-LAYER PWL PROFILES OF WIND SPEED 

The unperturbed state is assumed to be in hydrostatic 
equilibrium such that V(Pj + Pjgy) = 0, where Pj and 
Pj are the pressure and density, and j = a,w stand for 
air and water. The equilibrium velocity is approximated 
by a PWL profile, in which the liquid is at rest (so that 
the effect of the water drift velocity is neglected) : 

P = Pw,U = for y < 0, 
P = Pa-.U = Say for < ?; < L^, 
P^ Pa,U ^Ua for y > La, (1) 

where Ua is the wind speed at the reference height y — La 
in the air; La is the thicknesses of the shear layer in the 
air (that has the meaning of the roughness of the air- 
water interface); Sa is the wind shear. 

The stability properties of the equilibrium state (1) are 
investigated by considering the following normal modes: 



$(a;, y, t) — (j){y)exp{—iujt + ikx), 



(2) 



where $ stands for any of the perturbed variables and cj) 
is its amplitude; w = w^ -I- iuJi and C = Lo/k = Cr + iCi 
are the complex frequency and phase velocity, and fc is a 
real wave number. 

In the case of the PWL approximation for the veloc- 
ity in Eqs. (1), the dispersion relation rewritten in the 
present notations is as follows ([6|, [l6|): 

[uj'^iPw + Pa) - Pa^Sa - [pw ^ Pa)gk - (jk^] X 

X [C^ - SaF{kLa)] = -paUJ^Sa [GikLa) + 1] , (3) 

where 

F{x)]^ -exp{-2x)-- + x, G{x) =exp{-2x)-l, (4) 



F{x) ^x^ + 0{x^), Gix) = a; -^ 0{x^) for a; < 1, 



F{x) ^x + 0{x°), G{x) = -I + 0{exp) for x > 1, 

and 0{exp) denote exponentially small values in x 3> 1. 

Gaining a deeper understanding of the nature of the 

Miles and KH regimes of instability is facilitated by the 



introduction of the following two sets of characteristic 
scales (denoted by stars): 



(0 P* ^ Pw, 


L^ = La, i, = \/L^/g, 


(5) 


{a) p, = p^, 


L* = Lg, i, = y/Lt,/g. 


(6) 



where Lg = Ua/g is the gravitational length. As the 
characteristic properties of the two modes of instabil- 
ity will be unfolded in the next sections, it will become 
clearer that the Miles mechanism is better described by 
scaling the physical variables with the characteristic val- 
ues given in Eqs. (5), while the scales depicted in Eqs. 
(6) are naturally suited to describe the KH mechanism. 
For now, though, the physical variables are scaled in a 
general way as follows: 



k — kL^:, Lo — ujt^ 



9P 



■LV 



C = 



u = 






2 ^ Pa_ 



(7) 



Consequently, the dispersion relation (3) is given now 
in the following dimensionless form 



[u;''{l + ^)-e^QSa 



{l~e^)k 



at^ 



x[iD- SaF{kLa)] = -e^UJ^Sa[GikLa) + 1] , 



(8) 



and each of the particular cases (represented in Eqs. (5) 
or (6) is determined by fixing the values La or Ua, re- 
spectively, 



Ua 



:, Sa = Ua, W 



ii)La^l, Ua 

vy^a 

(a) La = -rr^, Ua = 1, Sa = ^— , W 



-'9' 



k = ka, (9) 

k = fcg.(lO) 



Thus, the system's stability is determined by the fol- 
lowing three dimensionless parameters: density ratio e^, 
surface tension a, and Froude number Fr — Ua/\/gLa 
{Ua = Fr and La = l/Fr"^ in Eqs. (9) and (10), respec- 
tively). 

Before turning to the solutions of the dispersion equa- 
tion (8) it is noted that the dimensionless surface tension 
coefficient may be expressed as a ratio of two character- 
istic lengths: 



Uc/F^, 



fll) 



where L^ — \/cr/{gPw) is the constant capillary length. 
Doing so it is evident that the dimensionless surface ten- 
sion coefficient is uniquely determined by the character- 
istic scale i* and vice versa. 

Dispersion relation (8) is a cubic equation for Qa as 
a function of ka, which may be solved analytically by 
employing the Cardano solution. However, such a so- 
lution results in complex expressions that conceal the 
physical meanings and tendencies of the various regimes 
of instability. Instead, in the following sections disper- 
sion relation (8) will be solved in various regimes of 



parameter space with the aid of an asymptotic expan- 
sion in the natural small parameter of the problem, i.e., 
e^ = Pal Pw ~ 10"'^. As will be seen, such expansions 
result in simple and physically transparent solutions. 



III. MILES' REGIME: SHORT WAVES 
GENERATED BY WEAK WINDS 

Dispersion relation in the leading order in e. — In or- 
der to investigate the Miles regime the dispersion relation 
(8) is investigated in terms of the characteristic scales in 
Eqs. (5) denoting the corresponding dimensionless values 
of parameters by subscript a. To simplify the resulting 
relations, it is convenient to introduce the following aux- 
iliary scaled variables: 



U = U;a/Ua, C=Ca/Ua. 



(12) 



Substituting definition (12) into the dispersion relation 
(8) yields 

Ua 

^£^[^3 + ^^ „ F{h){^^ -^ + m) + '^'G(fca)]. (13) 

First, the generation of short waves (with wave length 
A of the order of roughness La) is considered, such that 
all the relevant physical variables are of the order of e° 
with respect to the small density ratio 



fca~e°, t/a~e°, S^e°, 



(14) 



In that case, Eq. (13) has the following two solutions 
in leading order in e 



iO - F{ka) =0, w ^^- 



0. (15) 



To classify the waves that are described in Eq. (15), it is 
useful to return momentarily to the dimensional variables 
in Eqs. (15): 

to - SaF{kLa) = 0, p„w^ - pyjgk - ak^ = 0. (16) 

The second equation in (16) evidently describes the 
gravity-capillary waves within the current approximation 
of low density ratio, while the first equation in (16) rep- 
resents the inertial waves uj — kUa = kiC — Ua) = in 
the limit ka = kLa — > oo {F{k) « ka for ka — > oo). 
With this in mind, these waves will be called inertial and 
gravity-capillary waves, respectively. 

Both modes of wave propagation that are described 
in Eqs. (15) are stable to leading order. Furthermore, 
since all the coefficients in the dispersion relation (13) 
are real, the higher order corrections to Cu in Eqs. (15) 



reasoning breaks down, though, when the inertial wave 
merges with one of the gravity-capillary waves (which 
happens when two roots of Eqs. (15) are merged into 
one double root, see Appendix A). In that case, higher 
order corrections are of order of e",n = 1, 2..., and may 
lead to instability. The condition for a double root to 



occur is: 



Ua = 



\Jka -h gfcg 
F(ka) ' 



(17) 



where the resonance condition is resolved with respect to 
IJa- Thus the system can be driven to unstable state via a 
resonant interaction of the gravity-capillary wave on the 
water surface with the inertial air wave, i.e. the transi- 
tion from stable to unstable waves is associated with the 
merging of the two waves. The image of the resonance 
condition (17) in the Cardano solution of a cubic equa- 
tion is the vanishing in leading order of the discriminant 
of the cubic equation with real coefficients. As a result, in 
such a case small corrections of the order of e^ neglected 
in the solutions (15) can provide two complex conjunct 
roots and, hence, to destabilize the system. The real part 
of the double roots of Eq. (13) to leading order in e are 

^, =F(fc,), a=F(fca)/fca, (18) 

while the growth rate to leading order is (Appendix A) 



U^, ^ e\\\^F{ka)\G{ka) 



(19) 



The instability that occurs due to the resonance condi- 
tion (17) in flows with the two-layer PWL profile of wind 
speed is quite similar to Miles' resonance condition for 
the instability of flows with smooth velocity profiles (as is 
discussed in J16j). and the corresponding unstable mode 
is called Miles mode. The present asymptotic approach 
leads to the exact resonance condition (17) which pro- 
vides an explicit expression for the resonant wind speed 
that leads to instability. The resonant nature of the in- 
stability is confirmed by the numerical solution of the 
same dispersion relation by obtaining a smoothed near 
resonance solution in an asymptotically small e-vicinity 
of the resonant wave numbers 16| . 

Although the scaled values of the complex frequency 
UJ and phase velocity C, are independent of the surface 
tension, the unsealed complex frequency (Dq — ojUa and 
phase velocity Ca = CUa do depend on the surface ten- 
sion through Ua in Eq.(17), which strongly varies with 
a. The unsealed dimensionless frequency, phase velocity 
and growth rate are given by: 



= Jka + akl, 



Ca 



^^''U.^W^ii)±l. (20) 



2F{ka) 



are of orders e 



2ri 



1,2..., and are real. That line of 



The frequency, growth rate, phase velocity and the res- 
onant wind velocity are depicted for several values of 



a (or, equivalently, the characteristic roughness length 
i» = Lq, see relation (11)) as functions of the wave num- 
ber in Fig. [H 



TABLE I: System parameters at minimum wind velocity 
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FIG. 1: Dimensionless scaled: 

(a) frequency ujr = u!ar/Ua;{h) growth rate uji/e = ujat/i^Ua); 

(c) phase velocity Cr = Car/Ua; 

Dimensionless unsealed: 

(d)wind velocity tJa; (e) frequency Uar', (f) growth rate uJai/t; 

(g) phase velocity Car vs wave number ka ; the curves marked 

by 1, 2, 3, 4 correspond to ct = 0; 0.15; 0.5; 2.0. 

According to Fig. [1};, the resonant wind velocity de- 
pends on the wave number, or, physically more correctly, 
the resonant wave number depends on the wind velocity. 
For a nonzero surface tension {a > 0) Fig.[T]demonstrates 

a threshold (critical) wind speed Ua' for ka = ka' nec- 
essary to make the waves grow (see also Table 1). In 
particular, for any given supercritical wind and wave ve- 
locities the resonance condition is satisfied for two wave 
numbers (Fig.[TJl), which correspond to waves with differ- 
ent growth rates: longer waves (with smaller wave num- 
bers) grow faster, since the growth rate monotonically 
increases with the wave length (Fig. [If). Similarly the 

minimum gravity-capillary wave speed CaT is achieved 

at the wave number ka = ka that differs from that for 
the wind speed (see Figs. Id and Ig) 



ct'-y = n^, ^r^' = i 



(21) 



In addition, according to Fig.[l]the value of the roughness 
(or, equivalently, dimensionless surface tension) strongly 
influences the marginal stability values of the wind veloc- 
ity as well as the wave phase velocity, while the growth 
rate depends only slightly on the roughness (see Fig.[T]:). 
The critical values of the dimensionless as well as the 
corresponding dimensional parameters presented in Ta- 
ble I for several values of the dimensionless surface ten- 
sion a = 0.15; 0.5; 2.0 have been calculated at minimum 
wind velocity (ka = ka ) by using the following rela- 



{ka = ka ) for Miles' regime in air- water system. 
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tions: 

La = Lc/V<J, Ua = Ua\/gLa, Ca ~ Ca\/ gLa, 

ka^ka/La, \a^2^/ka, Lg = U^/g. (22) 

Analysis of results in Table L. — The critical value of 
the wave length A ~ (0.7 -^ 1.2) • lO^^m is of the order of 
the wind profile inhomogeneity (roughness), La '^ (0.2 -^ 
0.7) • lO^^m, as well as of the order of the gravity and 
capillary lengths Lg - (0.7 ^ 1.5) • 10"^^^ ^^ ^ o.27 • 



X 



Lr 



(23) 



The closeness of La and Lg follows from the following 
estimation for Ua 



(cr) 



Ua 



0.75-^1.6- 1 



Ua 



(24) 



The values of the critical parameter in Table 1 are close 
to those obtained numerically by [1^] for two-layer PWL 
profiles in Eq. (1). 

At this point it is natural to discuss the validity of the 
results obtained above for the PWL model for a smooth 
wind velocity profile. Since the wavelength is of the or- 
der of the roughness value in the critical region of pa- 
rameters A*-^''-' ~ La (Table 1), one should not expect a 
model with discontinuous profiles to be compared favor- 
ably with the smooth profiles being modeled, it rather 
models schematically the influence of the wind-shear ef- 
fects on the interface instability. Indeed, the growth rate 
coefficients in Fig. lb significantly overestimate the val- 
ues calculated numerically for smooth wind profiles in 
0. Nevertheless, the values of the wind velocity, rough- 
ness and wave length in Table 1 are in fair agreement by 
the order of magnitudes with experimental data for low 
wind velocities [l0| {La ~ 0.3 • 10~^m at Ua ~ 0.23?7i/s). 
The latter results are also consistent with the results of 
Miles' model for the smooth logarithmic profile in the 
limit of small density ratio: Ua ~ VgUc ~ 0.33?7i/s, 
Xi^r) ~ 1 7 . iQ-'^m. As was mentioned above, the two- 
layer PWL approximation leads to the present resonance 
model, which is quite similar to the Miles critical-layer 
resonance model for smooth velocity profiles [l6| . 

Region of asymptotic nonuniformity in small ka ^ £■ — 
The solution presented above loses its validity for suffi- 
ciently long waves Aq = X/La (small ka = 27r/Xa)- The 



nonuniformity of the asymptotic solution is exhibited in 
the infinite values of the resonant wind velocity Ua{ka), 
growth rate oJaiika), and phase velocity Carika) at low 
values of ka in Fig. 1. That resonant solution has been 
obtained by a formal asymptotic expansion in e that has 
been applied to the cubic dispersion equation (13), with 
the scaled eigenvalue w expanded as follows (Appendix): 



2,~, 



UJQ + €Ll!i + e bJ2 + 



(25) 



Such an expansion is a direct result of the assump- 
tions expressed in estimations (14), that both the wave 
number as well as the wind velocity are of zcroth leading 
order in e. However, as the wave lengths and the wind 
velocity increase, the asymptotic expansion (25) loses its 
validity. Thus, turning now to small values of the wave 
number ka ^ e (see Appendix A), and inspecting again 
dispersion relation (8), results in the following scalings of 
the remaining physical variables: 



Wo 



^e,^(n 



^ar 



,- . , /7 /"'' .- 


1 

i;3/2 

e 

i;3/2 


1 


1 1 /^ 


£3/2 
1 




V^' 



(26) 



Thus the region of nonuniformity indeed corresponds to 
the longer waves with ka ^ £ and higher wind velocity 
tJa ^ e~3/2 fj[^g^j^ [^ Tfjg^ assumed in estimations (14). 
To further develop the asymptotic solution in the re- 
gion of nonuniformity, it is recalled that the classical KH 
mode (vortex sheet instability) has the following orders 
in terms of the characteristic length Lg ( e.g. 14|): 




e, Ug ^ 1. 



(27) 



Rewriting relations (26) in terms of the characteristic 
length Lg yields by using estimations (26)-(27): 




Comparing Eqs. (27) and (28) results in the following es- 
timate of the characteristic roughness for which the res- 
onant solution describes the classical KH mode: 



La — La/ Lg ~ e 



(29) 



This is also consistent with the estimation (26) for Ua — 
Ua/\/gLa ~ e~3/2 (accountiug for the definition of the 
scaled dimensionless roughness length La = Lag/U^)- 
Substituting now (29) into (28) leads indeed to the 
asymptotic scales already defined above in relations (26). 



In particular, according to estimates (26) the growth rate 
of those moderately-short waves is expected to be much 
larger than that (see Eqs. (20)) for the Miles short waves. 
The term moderately-short-wave mode stems from the 
fact that its wave length is much larger than that cor- 
responding to the scales of the interface roughness and, 
simultaneously, much less than the scale of gravity waves: 

La < A < Lg 

or in dimensionless form using estimations (28)- (29): 



1 X X 1 Lg 1 

La e La e-^ 



(30) 



Finally note that the assumption (14) for the asymp- 
totic orders of the values correspond to the outer solution 
in terms of the method of matched asymptotic expan- 
sions [20] , while within the region of nonuniformity (26) 
or (27) of the corresponding outer asymptotic solution 
(Miles' mode) the inner solution should be constructed. 



IV. KH REGIME: MODERATELY-SHORT 
WAVES GENERATED BY STRONG WINDS 

Dispersion relation in the leading order in e. — The 
asymptotic solution is now constructed in the inner re- 
gion of small ka ^ e, where the outer solution described 
by relations (17) and (20), is not uniformly valid. Follow- 
ing estimates (26) it is useful to introduce scaled variables 
of zeroth order in the inner region: 



e' u;a 



Ca=CaV~^--e", Ua 



= Uae'/' 



(31) 



Substituting the inner variables (31) into dispersion 
relation (8), and expanding F{ka) and G{ka) in Taylor 
series in A:a ^ 1 by using (4), yields in the leading order 
in e the following reduced dispersion relation: 

ul - Uakl^l + UJaKiUaka - 1) + Uakl = 0. (32) 

In order to estimate the effect of the surface tension the 
experimentally observed roughness value La ~ 10"'^ is 
adopted for 25to/s < Ua < 50m/s (Fig. 3b in 21]). 
Hence, since a = L^/Ll ~ 10 {Lc ~ 3 • 10~^), the sur- 
face tension term in the dispersion relation (32) that is 
~ e^(T has been neglected in comparison with the remain- 
ing terms ^ e°, and hence does not appear in Eq. (32). 
Consequently, accounting for the surface tension effect in 
the classical criterion of the KHI has no physical mean- 
ing. 

Gravitational length as the characteristic scale. — For 
comparison convenience with the limit of vortex sheet 
instability (i.e., the classical KHI limit) the dispersion 
relation (32) is rewritten in terms of the characteristic 
scales in Eqs. (6) based on the gravitational length Lg = 
Ua/g and scaled by e: 



OJn 



Lak LU + U)gkg{kg - 1) + Lat = 0, 



(33) 



TABLE II: Typical parameters of the air-water system at 
hurricane winds for KH mode. 



curve La ~ La (kg) or the value La = 0, but below the 
marginal stability curve La = La {kg): 



Ua 


kg A - 


= 27re2 


9kg 


La 


La 


_ La a 
7^ 


m/s 




m 




m 






25 


0.2 


2.5 




0.001 




0.4 


25 


1.0 


0.5 




0.001 




0.4 


25 


5.0 


0.1 




0.001 




0.4 


50 


0.2 


10. 




0.001 




0.1 


50 


1.0 


2.0 




0.001 




0.1 


50 


5.0 


0.4 




0.001 




0.1 



;(i) 



•-(2 



or, equivalently, in terms of a complex phase velocity 



(73 - LakgC^, + Cg{l - fc-1) + La = 0, (34) 



where 



^9 = ^ C'ff 



e" iOa 



e'^LUn 



La — e La 



(35) 



For such a choice of the characteristic scales, the disper- 
sion relation (33) or (34) explicitly contains the dimen- 
sionless roughness La — La/Lg, so that the latter may 
be set to zero, and thus obtain the dispersion relation for 
the classical KHI limit (see Eq. (39) below). 

In Table II the wave lengths A are estimated at wind 
speeds Ua = 25 and 50m/ s {La ~ IQ-^m, [HI) for the 
wave numbers kg = 27re^C/^/(gA) of the order of e°. It is 
seen that the typical wave lengths in the KH regime are 
ranged from A ~ 0.1m at Ua ~ 25m/s to A ^ IOttt. at 
Ua — 50m/s. 

General relations for marginal stability curves. — The 
cubic equation (34) with real coefficients has either one 
real and two complex conjugate roots, or three real roots 
at least two of which are equal, or three different real 
roots. The various cases depend on the value of the dis- 
criminant Q that may be either positive or zero, or nega- 
tive, respectively. It is convenient to write the condition 
of zero discriminant in the following way: 



aLa'^ + [3La^ +7 = 0, 



(36) 



where 



a{kg) = 4fc3, P(kg) - kl-20kg^^, j{kg) - -4(1-^^1)3. 

Consequently, the corresponding (physically meaningful, 
i.e. only positive) marginal roughness lengths are 



L^:''\kg) = 



-/? ± V/?^ - 4a7 



2a 



(iW<i(2)). (37) 



In general, Eq. (36) has either two or one positive so- 
lutions (37). The number of solutions depends on the 
existence of a real solution for La , while the instabil- 
ity occurs in one of the following two cases: the re- 
gion of instability lies above either the marginal stability 



L'a'' < La <L'a'> OyO< La <L'^a>- 

Consider now the marginal stability frequency ujgm 
that is obtained by setting the imaginary value of cj to 
zero {ujgi = 0) . This may be equivalently obtained by set- 
ting the discriminant Q = equal to zero. The marginal 
frequency and phase velocity are then given by: 



Wo 



where 



9 9^^' ^ gyn 



q^^j_Lakg_ 

2 ^ 27 






kg -I 



(38) 
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The dimensionless values of the characteristic rough- 
ness, frequency, growth rate, real and imaginary phase 
velocities vs the wave number are depicted in Figure 2. 
The marginal curves 1 and 2 set the boundaries of the 
instability region in Figs. 2a-2e. It will be shown briefly 
in the next subsection that the marginal stability curves 
1 and 2 represent KH-like and inertial-like modes, re- 
spectively. In Figure 2a the marginal stability curve 1, 
La — La {kg), corresponds to the classical KH mode gen- 
eralized by accounting for the wind-shear effects. The ap- 
pearance of the second marginal curve 2, La = La {kg), 
which is absent in the classical KH mode, is entirely 
due to the wind-shear effects. At constant values of the 
roughness its influence, as seen from Figs. 2a and 2c, 
leads to cut off of the growth rates in the short wave 
limit. Curves 3 and 4 depict the values of the complex 
frequency and phase velocity for La = 0.5 and La = 0.12, 
which correspond to La ~ 10~^m, {Ua = 25m/s and 
Ua = 50m/ s, [21|). In Figures 2c and 2e the growth- 
rate curves 5 (imaginary frequency and phase velocity) 
are depicted for the classical KH approximation for the 
same wind speeds Ua = 25m/ s and Ua = 50m/s. It 
is seen that they are very close to those for the PWL 
approximation in the whole range of the wave numbers 
at Ua = 50m/ s, but only at slightly supercritical wave 
numbers at Ua = 25m,/ s. 

KH and inertial limits for marginal curves. — The limit 
La ^ corresponds to the classical KH mode (see Fig. 
2a), and the point fc^ = 1 where the marginal curve 1 
crosses the fc^-axis is the threshold of the classical KHI, 
so that the instability occurs for kg > 1. Indeed, the 
dispersion relation (33) may be additionally expanded in 
I/a ^ 1, and thus reduced to the limit of the classical 
KH mode (vortex sheet instability) : 



+ U}gkg{kg - l) = 0, 



(39) 



where the trivial solution Ug — should be omitted. To 
make the comparison evident the dispersion relation (39) 
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FIG. 2: Marginal stability and growth rate curves vs wave 
number, (a) roughness La\ (b) frequency Ugr', (c) growth rate 
ujgi\ (d) real phase velocity Cgr', (e) imaginary phase velocity 
Cgi. the KH-like and inertial-like marginal curves 1 and 2 
correspond to signs - and + in Eq. (37); the straight line 2' 
is the asymptote Cgr = 1 at fcg — » oo. Dashed curves 3 and 
pointed curves 4 depict growth rates for La ~ 0.5, La = 0.12 



Ua, equality (43) may be rewritten in dimensional form 
'AS C — Ua- Hence, the corresponding marginal curves 2 
in Fig. 2 are termed inertial-like. 

Definition of the KH-like regime of instability. — Fol- 
lowing the discussion in the previous sections, the region 
between the marginal curves 1 and 2 in Fig. 2, is called 
the KH-like regime or finite-shear-KH regime. Such a def- 
inition seems to be physically plausible and convenient 
in order to distinguish the KH-like mode that is char- 
acterized by moderately-short waves and strong winds 
from Miles' mode that is characterized by short waves 
and weak winds, as well as by the resonant nature of the 
instability. 

Note that the high-fc^ limit of the upper marginal curve 
2 in Fig. 2d in the region of the KH-like instability is 
matched asymptotically with the axis Cg = for the 
Miles' mode at e —+ 0. To see that, according to the 
conventional asymptotic matching procedure of the inner 
and outer expansions [20], the inner (KH-like) marginal 
stability curve Cg = 1 in (43) should be rewritten in the 
outer (Miles) limit. Then by using relations (28) and (35) 
between the inner and outer variables kg , Cg and ka, Ca, 
one has that the region of the KH-like instability shrinks 
in the outer Miles limit since the marginal curve 2 in Fig. 



at Ua = 25m/s, Ua = 50m/s respectively {L = lO" m). ^ approaches zero in that hmit (fc, 



Point-dashed curves 5 in Figs. 2c and 2e depict growth rates 
for the classical KH approximation ujgi = -J k^ — kg 



oo, e 



0)- 



Comparison with the numerical results in Iq]. — The 
results presented in Fig 2a have been reconstructed in 



the plane (, 



kLa) for comparison with the numcr- 



may be rewritten in dimensional form as: 

Pw 



(40) 



ical results presented in [y|. The marginal curves ob- 
tained for the present asymptotic solutions in e in Fig. 
3, are very close to those calculated numerically (see Fig. 
3 in |6[). According to (41) the upper marginal curve 



Thus, a vicinity of the fcg-axis in Fig. 2a is the region of 
applicability of the classical KH approximation, in par- 
ticular, the region includes a part of the marginal curve 
1, where La <^ 1, and the marginal curves 1 are termed 
KH-hkc. 

On the other hand, the condition La ^ 1 is only neces- 
sary but not sufficient for the convergence of the KH-like 
theory to the classical KHI limit. Indeed, according to 
Eqs. (36) in the limit of fcg ^ oo the roughness length is 
given by: 



f(2) 



vanishes as 2 /kg for kg 



•-{2) 



OO {kgLa —> 2). Using 



relations (28) and (35) this yields the constant asymp- 
totic limit for the upper marginal curve kLa — 2e for 
kLg = oo. The point of touching the bottom marginal 



curve in Fig. 3 



e , kLa = 0) corresponds to 



ii^) « 2/kg 



0, 



while the dispersion relation (34) yields: 



2kgUjl + LJgk^ 



0. 



(41) 



(42) 



Equation (42) is indeed quite different from Eq. (39) 
for the classical KH mode. The solution of Eq. (42) in 
the limit of fcg ^ oo (represented by the dashed straight 
line 2', in Fig. 2d) is given by: 



Co 



1. 



(43) 



Remembering now that the phase velocity has been 
made dimensionless by using the characteristic velocity 



the point {kg = 1, Lq = 0) on the bottom of the marginal 
curve in Fig. 2a. Thus, having in mind that the region 
characterized by low values of La , and fcg > 1 is the re- 
gion of applicability of the classical KH mode, it may 
be concluded that the point is the critical point for the 
classical KH mode rather than the point of 'a continuous 
transition from the Miles-type instability to the KH-type 
instability' (cited from [6|). 

The closeness of the marginal stability curves obtained 
numerically by [6[ for arbitrary La to those in the present 
study, in spite of the fact that the latter has been selected 
on the basis of the asymptotic scaling hypothesis (29) for 
the classical KH regime for small La, supports the result 
dispersion relation (33) 



V. CONCLUSIONS AND DISCUSSION. 

Asymptotic expansion in the small air-to-water mass 
densities ratio e^ <C 1 that has been applied to the two- 



^ 




FIG. 3: KH-like and inertial-like marginal curves 1 and 2, kLa 
vs \JkLg (for comparison with the numerical results [Q] for 
zero water drift speed). The instability region is between the 
marginal curves 1 and 2. The asymptotic abscissa \Jk,Lg — 
e~^ of the touching point kLa = of the bottom marginal 
curve is the threshold point of the classical KHI (kg = \,La ~ 
0); kLa = 2e is the asymptotic limit for the upper marginal 
curve for kLa ~* oo. 



layer PWL approximation of wind profile (the simplest 
modeling of wind-shear elTects) provides an explicit and 
physically transparent solution of the dispersion relation. 
Based on the conventional method of matched asymp- 
totic expansions 2Q\, it has been shown that the Miles 
and KH regimes are described by the outer and inner 
asymptotic solutions in small e, respectively, and operate 
at quite different scales of the unperturbed wind speeds 
and perturbation wave lengths: 

(i) The Miles mechanism is relevant for weak winds, 
which generate short waves with wavelengths A of the 
order of the scales of the wind profile inhomogeneity 
(roughness) La, as well as the characteristic capillary and 
gravitational length, Lc and Lg 



\ 


Lc 


Ln 










9 n 


^f^ 


Ft ^ 


La 


La 


La 






La 





JLa 




a 




" e , ^i 


\ 


~ (-, 





9 




y 9 




^JLag 




6.(44) 



(ii) The KH mechanism is relevant for strong winds 
which generate moderately-short waves with wavelengths 
A much larger than the roughness La and the capillary 
length Lc, but much shorter than the gravitational length 



Lg {La ^ Lc <t: X <$: Lg 




According to estimates (44) and (45) for Froude num- 
bers the characteristic wind speed in the KH-like regime, 
Ua ^ e~^^'^y/gLa, is much larger than that, Ua "^ 
e'^y/gLa, in the Miles regime. In addition to the large 
difference in the scales of the perturbation wavelength 
and wind velocity for the two modes, the growth rate 
of the KH mode is much larger than that for Miles' 
mode. Furthermore, accounting for the surface tension 
effect for prediction of the minimum wind speed is cor- 
rect for Miles' regime, and has no physical meaning in 
the classical KH model. Estimates (44) and (45) for the 
characteristic wave lengths in the Miles and KH regimes 
are consistent with respect to the orders of magnitude 
with those in Table I and Table II, respectively. In par- 
ticular, the typical wave lengths in Miles' regime vary 
from A ~ 0.2m for Ua = 0.26cm/s to A ~ 0.7cm for 
Ua — 0.58cm/ s, while in the KH regime they range from 
A ^ 0.1m for Ua — 25m/s to A ~ 10m for Ua = 50m/s. 

An important question is to what extent the PWL ap- 
proximation provides a realistic representation of smooth 
wind profiles for both the Miles as well as the KH regimes. 

Miles' regime. — A reasonable similarity with smooth 
wind profiles with regard to the stability characteristics 
should not be expected for the A4iles modes in the two- 
layer PWL approximation of the wind velocity with wave 
lengths of the order the roughness scale, it rather mod- 
els schematically the influence of the wind-shear effects 
on the interface instability. However, even in such an 
approximation, the boundaries of the instability region 
in the Miles' regime are in fair quantitative agreement 
with more exact theories, and only the growth rate value 
is significantly overestimated. For small density ratio 
the specific feature of the PWL approximation in the 
Miles regime is the resonant nature of the solution that 
is characterized by a spectrum with a finite number of 
localized almost-discrete wave lengths. Indeed, numer- 
ical solutions of the two-layer PWL approximation for 
the wind speed demonstrate that the instability occurs 
in narrow intervals that are localized about either one 
(for CT = 0) or two (for a > 0) resonant wave numbers, 
which are determined explicitly in the present modeling 
in the limit of small e. Furthermore, the number of reso- 
nant unstable waves grows with the number of the PWL 
sections N such that the discrete spectrum tends to the 
continuous spectrum as the PWL profile tends to the con- 
tinuous one in the limit of large N. The convergence rate 
strongly depends on the system parameters. In partic- 
ular, numerical solutions of the Rayleigh equation with 
A^-layer PWL approximation for the wind profiles [91], in 
the Miles regime, indicate that the convergence in N is 
very slow for small density ratios e^ « 10^'^, but quite 
good for larger values e^ ~ 10~^. 

KH-like regim,e. — Since the details of the wind profile 
are less essential when the wavelengths are much larger 
than the profile inhomogeneity length scale, it may be 
naturally expected that the two-layer PWL approxima- 
tion of the smooth profiles should be better within the 
moderately-short KH regime than for the short-waves 
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FIG. 4; Growth rate vs wave number in the plane (a) {ka,iOai) 
and (b) (kg,LJgi); curves 1, 1' are for Ua = SOtti/s, curves 2, 
2' are for Ua = 2bm/s {La = 10"''). Curves 1 and 2 are for 
tanh-profile of the wind curves 1' and 2' are for the two-layer 
PWL approximation of the wind profile. 



Miles regime. Numerical results for the smooth profile 
U{y) = tanh{y) have been obtained by solutions of the 
well-known Rayleigh equation for the inviscid flows (e.g. 



0). Indeed, unlike the behavior within the Miles regime, 
the two-layer PWL approximation of the wind velocity 
for the KH regime provides a fair quantitative agreement 
for the growth rates with those for the smooth wind pro- 
file (Fig. 4). 

Figure 4 depicts the growth rate for the same parame- 
ters as Fig. 2c. As it is seen in Fig. 2c, the classical KH 
approximation of the growth rates is in a fair quantita- 
tive agreement with those in the two-layer PWL approx- 
imation (and hence with the results for U = tanh{y)) 
for sufficiently large wind speeds. The growth rates de- 
picted in Fig. 4a in the characteristic scales in Eqs. 
(5) typical for Miles' region of parameters have been 
reconstructed in Fig. 4b in the characteristic scales 
in Eqs. (6) natural for the KH regime. Note that if 
kai < ka2 it may be that fcgi > kg2 if kiUi > fc2C^2 since 
kgi — e^kxUl/g > fcg2 — e^feC/|/g as it occurs for the 
locations of the maximum growth-rates in Fig. 4. 
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APPENDIX A. ASYMPTOTIC RELATIONS FOR 
MILES' REGIME 



Substituting now expansion (A4) into (Al) yields the fol- 
lowing new first two terms in expansion (A2): 



Omitting for simplicity subscript a at all parameters 
and bar under k, U and a in dispersion relation (13) 
yields 



$(w) = -e2^(w), 



(Al) 



where 



^(Cb) 



^(J}) 



F{k)] 
k 



(2.2 



k + aP 
t/2 



)^u^ + —u- F{k) (^^^ + — -Co)+ Co^G{k) 



Equation (Al) leads to the following formal asymptotic 
expansion 



(I) = a)o + £^^^2 + 



(A2) 



Substituting (A2) into (Al) yields for the first two terms 
in expansion (A2): 



$(a)o) =0, (1)2 = - 



^(wo) 



9$(aJo)/5e 



2' 



(A3) 



The formal solution (A3) is vahd if d^{uJo)/de^ ^ 0. 
Otherwise (the resonance case) the asymptotic expansion 
(A2) should be replaced by: 



u) — Cji^ + ecDi + e L02 + ... 



(A4) 



Wi 



2*(c^o) 



92$(c;)o)/ae2 



^(cZ-o) = 0, 



*(^o) = F^{k)[l + G{k% a2$(Jo)/5e2 =. 4F{k). 

As a result, the unsealed value lui = loiU with U given 
by Eq. (17), is equal to 



, 1 + G(fc),, 



(A6) 



The asymptotic expansion (A4) fails for small fc ~ e, 
when the leading order term in Eq. (A4), Cjq = F{k) ~ 
e°fc2 becomes of the same order as the next term euji ^ ek 
(here the estimates (4) for F{k) and G{k) have been 
used). Returning to the unsealed variables (12) and en- 
tire notations, this yields for the region of nonuniformity 
of Miles' mode: 



e, Wo 



Vk ^ y/e, 



euJi 



(A7) 



